Electronic correlations at the a— 7 structural phase transition in paramagnetic iron 
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We compute the equilibrium crystal structure and phase stability of iron at the a(bcc)-7(fcc) 
phase transition as a function of temperature, by employing a combination of ab initio methods 
for calculating electronic band structures and dynamical mean-field theory. The magnetic correla- 
tion energy is found to be an essential driving force behind the Q-7 structural phase transition in 
paramagnetic iron. 
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The properties of iron have fascinated mankind for 
several thousand years already. Indeed, iron has been 
an exceptionally important material for the development 
of modern civilization and its technologies. Neverthe- 
less, even today many properties of iron, e.g., at high 
pressures and temperatures, are still not sufficiently un- 
derstood. Therefore iron remains at the focus of active 
research. 

At low pressures and temperatures iron crystallizes in a 
body-centered cubic (bcc) structure, referred to as a-iron 
or ferrite; see Fig. [TJ In particular, at ambient pressure 
iron is ferromagnetic, with an anomalously high Curie 
temperature of Tc ^ 1043 K. Upon heating, iron ex- 
hibits several structural phase transformations 0, Q : at 
~ 1185 K to the face-centered cubic (fee) phase (7-iron 
or austenite), and at ~ 1670 K again to a bcc structure 
((5-iron). At high pressure iron becomes paramagnetic 
with a hexagonal close packed structure (e-iron). 

Density functional theory (DFT) in the local spin den- 
sity approximation gives a quantitatively accurate de- 
scription of the ordered magnetic moment and the spin 
stiffness of bcc-Fe Q , but predicts the nonmagnetic fee 
structure to be more stable than the observed ferromag- 
netic bcc phase 0] . Only if the spinpolarizcd generalized- 
gradient approximation (GGA) [5| is applied does one 
obtain the correct ground state properties of iron . 

Stoner theory of ferromagnetism 0] can give a qualita- 
tively correct description of several magnetic and struc- 
tural properties of iron, but predicts a simultaneous mag- 
netic and structural change at the bcc-fcc phase transi- 
tion with a local moment collapse while, in fact, the bcc- 
fcc phase transition occurs ^ 200 K above Tc; see Fig. [1] 
Clearly, to account for finite temperature effects of itin- 
erant magnets one requires a formalism which takes into 
account the existence of local moments above Tc- While 
the spin-fluctuation theory, which describes the param- 
agnetic state above Tc as a collection of disordered mo- 
ments, gives an overall good qualitative explanation of 
the pressure-temperature phase diagram of iron Q it fails 
to provide a reasonably quantitative description and, in 
particular, predicts the bcc-fcc phase transition to occur 
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FIG. 1: (color online) Schematic temperature-pressure phase 
diagram of iron (see text). 



below Tc- 

The LDA-hDMFT computational scheme Q, a com- 
bination of the DFT in the local density approximation 
(LDA) with dynamical mean-field theory (DMFT) [13], 
goes beyond the approaches discussed above since it ex- 
plicitly includes many-body effects in a non-perturbative 
and thermodynamically consistent way. LDA-I-DMFT 
was already used to calculate the magnetization and the 
susceptibility of a-iron as a function of the reduced tem- 
perature T/Tc Uj. The calculations gave overall good 
agreement with experimental data. The problem has 
been recently revisited by Katanin et al. |12l | who found 



that the formation of local moments in paramagnetic a- 
Fe is governed by the Cg electrons and is accompanied 
by non-Fermi liquid behavior. This supports the results 
obtained with the s-d model for the a-phase of iron [31 . 
A recent implementation of the LDA/GGA-I-DMFT 



scheme in plane-wave pscudopotentials [3, lE] now al- 
lows one to investigate correlation induced lattice trans- 
formations such as the cooperative Jahn- Teller distortion 
in KCuFa and LaMnOs. The method was not yet used to 
study structural phase transitions in a paramagnetic cor- 
related electron system with temperature (or pressure) 
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involving a change of symmetry. This will be the goal of 
the present investigation. 

In this Letter we employ the above-mentioned imple- 
mentation of the LDA/GGA+DMFT scheme [11 to 
explore the structural and magnetic properties of para- 
magnetic iron at finite temperatures. In particular, we 
will study the origin of the a-7 structural phase trans- 
formation, and the importance of electronic correlations 
for this transition. We first cornpute the nonmagnetic 
GGA electronic structure of iron [16[ . To model the bcc- 
fcc phase transition we employ the Bain transformation 
path which is described by a single structural parame- 
ter c/a, the uniaxial deformation along [001] axis, with 
c/a = 1 for the bcc and c/a = \f2 for the fee struc- 
ture. Here the lattice volume is kept at the experimental 
volume of a-iron (a = 2.91 A) [3 in the vicinity of the 
bcc-fcc phase transition, while the c/a ratio is changed 
from 0.8 to 1.6. Overall, the GGA results qualitatively 
agree with previous band-structure calculations [^. In 
particular, the nonmagnetic GGA yields the fee struc- 
ture to be more energetically favorable than the bcc one 
(see Fig.©. 

Next we apply the GGA-fDMFT approach [H, [l5| 
to determine the structural phase stability of iron. For 
the partially filled Fe sd orbitals we construct a basis 
of atomic-centered symmetry-constrained Wannier func- 
tions [TB]. The corresponding first-principles multiband 
Hubbard Hamiltonian has the form 
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where fii. 



creates (de- 
stroys) an electron with spin a in the Wannier orbital 
m at site i. Here Hgga is the effective low-energy 
Hamiltonian in the basis of Fe sd Wannier orbitals. The 
second term on the right-hand side of Eq. [1] describes 
the Coulomb interaction between Fe 3d electrons in the 
density-density approximation. It is expressed in terms 
of the average Coulomb repulsion U and Hund's rule ex- 
change J. In this calculation we use U = 1.8 eV which 
is within the theoretical and experimental estimations 
- 1 - 2 eV and J = 0.9 eV [TtI. Further, Hue is a dou- 
ble counting correction which accounts for the electronic 
interactions already described by the GGA (see below). 

In order to identify correlation induced structural 
transformations, we calculate [3] the total energy as 
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where Egga[p\ denotes the total energy obtained by 
GGA. Here (Hgga) is evaluated as the thermal aver- 
age of the GGA Wannier Hamiltonian. The third term 



on the right-hand side of Eq. ([2]) is the sum of the Fe 
sd valence-state eigenvalues. The interaction energy, the 
4-th term on the right-hand side of Eq. (2), {fiimanim'a') 
which is calculated in DMFT. The double-counting cor- 
rection Edc = 5E«mm',<T<T' C^mm'("«mcT)(n™v') Corre- 
sponds to the average Coulomb repulsion between elec- 
trons in the Fe 3c? Wannier orbitals calculated from the 
self-consistently determined local occupancies [T^]. 

To solve the realistic many-body Hamiltonian (1) 
within DMFT we employ quantum Monte Carlo (QMC) 
simulations with the Hirsch-Fye algorithm [l^ . The cal- 
culations for iron are performed along the Bain transfor- 
mation path as a function of the reduced temperature 
T/Tc- Here Tc corresponds to the temperature where 
the spin polarization in the self-consistent GGA-I-DMFT 
solution vanishes. We obtain Tc ^ 1600 K which, given 
the local nature of the DMFT approach, is in reasonable 
agreement with the experimental value of 1043 K and 



also with earlier LDA-I-DMFT calculations [11| . We find 
that Tq depends sensitively on the lattice distortion c/a. 
It has a maximum value for the bcc (c/a = 1) structure 
and decreases rapidly for other values. In particular, for 
all temperatures considered here the fee phase remains 
paramagnetic. 

In Fig. [2] we show the variation of the total energy 
of paramagnetic iron with temperature along the bcc-fcc 
Bain transformation path. The result exhibits two well- 
defined energy minima at c/a = 1 (at low temperature) 
and c/a = \pl (at high temperature), corresponding to 
the bcc and fee structures, respectively. We find that 
for decreasing temperature the inclusion of the electronic 
correlations among the partially filled Fe "id states con- 
siderably reduces the total energy difference between the 
a and 7 phases. In particular, the bcc-to-fcc structural 
phase transition is found to take place at Tgtruct ~ 1.3 Tc, 
i.e., well above Tc [1^. Our result for AT = Tgtruct — Tc, 
the difference between the temperatures at which the 
magnetic transition and the structural phase transition 
occur, is in remarkable agreement with the experimental 
result of AT ^ 200 K. This finding differs from conven- 
tional band-structure calculations which predict the mag- 
netic and structural phase transition to occur simultane- 
ously. Both Tgtruct and Tc vary sensitively with the value 
of the Coulomb repulsion U employed in GGA+DMFT 
calculation. We find that Tgtruct increases for increasing 
U values, whereas Tc decreases, in agreement with the 



Kugel-Khomskii theory [21 



In addition, we performed LDA+DMFT calculations 
to determine the phase stability of iron at the bcc-fcc 
phase transition as a function of temperature. In contrast 
to the standard band structure approach where it is es- 
sential that the spin-polarized GGA is used to obtain the 
correct ground state properties of iron, we find that both 
the LDA+DMFT and GGA-hDMFT schemes give qual- 
itatively similar results. In particular, both schemes find 
the bcc-to-fcc structural phase transition at ~ 1.3 Tc, 
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TABLE I: Calculated lattice constant a for the bcc lattice, 
volume V and bulk modulus B for the equilibrium phase of 
paramagnetic iron as a function of T/Tc- 



T/Tc 


Eq. phase 


a, A 


V, &v? 


B, Mbar 


(GGA) 


bcc 


2.757 


70.71 


2.66 




fee 


2.737 


69.20 


2.82 


0.9 


bcc 
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80.64 


1.48 


1.2 


bcc 


2.883 
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1.50 


1.4 


fee 


2.861 


79.03 


1.61 


1.8 


fee 


2.862 


79.13 


1.59 


Exp[l, 2, 22] 


bcc/fcc 


2.88-2.91 




1.62-1.76 
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FIG. 2: (color online) Variation of the total energy of param- 
agnetic iron obtained by GGA and GGA-hDMFT(QMC) for 
different temperatures. The total energy is calculated along 
the bcc-fcc Bain transformation path with constant volume 
(a = 2.91 A for the bcc phase). Error bars indicate the sta- 
tistical error of the DMFT(QMC) calculations. 

i.e., well above the magnetic transition. Explanations of 
the bcc-fcc structural phase transition and the fact that 
^struct 7^ Tc obviously need to go beyond conventional 
band structure theories. This clearly demonstrates the 
crucial importance of the electronic correlations among 
the partially filled Fe 2>d states. 

Next we perform a structural optimization and com- 
pute the equilibrium volume and the corresponding bulk 
modulus of paramagnetic iron (see Table H]) . The bulk 
modulus is calculated as the derivative of the total energy 
as a function of volume. We find that at the bcc-fcc phase 
transition the equilibrium lattice volume simultaneously 
shrinks by ~ 2 %, a result which is in good agreement 
with the experimental value of ^ 1 % [l|. The volume 
reduction is accompanied by an increase of the calculated 
bulk modulus. Overall, the equilibrium volume and bulk 
modulus computed by GGA-I-DMFT agree well with the 
experimental data [l|, d, . 

Finally we compute the square of the instantaneous 
local moment {ml) = ((Z]m["-™t ~ ^mil)^) of param- 
agnetic iron for the distortions c/a considered here. In 
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FIG. 3: (color online) Variation of the square of the local mag- 
netic moment calculated by GGA-f DMFT for paramagnetic 
iron. 

Fig. 13] wc show the result plotted for various tempera- 
tures. At low temperatures, the squared local moment 
depends quite strongly on the value of c/a, and is max- 
imal in the bcc and minimal in the fee phase, respec- 
tively. As expected, above T^, the square local moment 
gradually increases with temperature and becomes essen- 
tially independent of c/a, as indicated by the curve for 
T = 3.6 Tc in Fig. [3] (we note that this is only a hy- 
pothetical curve since at such an elevated temperatures 
iron is already in its liquid state). This finding has im- 
portant implications for our understanding of the actual 
driving force behind the bcc-to-fcc paramagnetic phase 
transition. For this we note that the squared local mo- 
ment (m^) determines the magnetic correlation energy 
— jl{ml), which is an essential part of the total cor- 
relation energy of the Hamiltonian (1). At high temper- 
atures, when the local moment is almost independent of 
c/a and the GGA-I-DMFT approach finds the fee phase 
to be stable, the contribution of the magnetic correlation 
energy to the bcc-fcc total energy difference is seen to 
be negligible. This changes markedly when the tempera- 
ture is lowered. Namely, upon cooling the contribution of 
the magnetic correlation energy gradually increases and 
becomes strong enough to overcome the DMFT kinetic 

energy loss Ek^n = Egga[p] + (Hgga) - T^m.k '^?r^k' 
the bcc phase as compared with the fee phase. Thereby 
the bcc phase with its larger value of the local moment 
is stabilized at T < 1.3 Tc- We therefore conclude that 
the bcc-to-fcc paramagnetic phase transition is driven by 
the magnetic correlation energy. 

In conclusion, we employed the GGA-fDMFT many- 
body approach to compute the equilibrium crystal struc- 
ture and phase stability of iron at the bcc-fcc transition. 
In particular, we found that the bcc-to-fcc structural 
phase transition occurs well above the magnetic tran- 
sition, and that the magnetic correlation energy is essen- 
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tial to explain this structural transition in paramagnetic 
iron. The above result and those for the equilibrium lat- 
tice constant and the variation of the unit cell volume at 
the bcc-fcc phase transition agree well with experiment. 
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